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Abstract 



We have searched for the decay B — > Z)^(2536)X and measured an upper 
limit for the inclusive branching fraction of B(B — > .CT^X) < 0.95% at the 
90% confidence level. This limit is small compared with the total expected 
B — > D^D^KX rate. Assuming factorization, the decay constant is 
constrained to be / D + < 114 MeV at the 90% confidence level, at least 2.5 
times smaller than that of Df. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 



One of the outstanding issues in B meson physics is the semileptonic branching fraction 
puzzle. Experimentally B(B —>■ Xiv) is measured to be (10.43 ± 0.24)% []l|, whereas theo- 
retical calculations have difficulties accommodating a branching fraction below ~ 12.5% [0. 
One way to reduce the theoretical expectations is through a two-fold enhancement in the 
assumed b —>■ ccs rate ||, which is estimated to be ~ 15% from the measured inclusive rates 
for B -> D+X and B -> ipX. 

Recently, Buchalla et al. Q and Blok et al. || have suggested that a significant fraction 
of the b — > ccs transition hadronizes into B — > DDKX. This is supported by CLEO's 
observation of "wrong-sign" Z) mesons from 5 decays, — > DX) = (7.9 ± 2.2)%, where 
the D comes from the virtual W + — > cs. The ALEPH |7j and DELPHI || collaborations have 
also observed sizeable B — > D^*'D^*'X decay rates. Exclusive .B decays involving wrong- 
sign D mesons can result from (1) resonant B — > D^*'D** decays, where the W + — > cs 
hadronizes to an excited Z)+ meson that decays into DKX; and (2) non-resonant I? — > 
£?(*)£) (*)/£" decays. This paper explores one possibility in the first case, namely, the decays 
i? — > Z) S ' 1 (2536)X where Z)^ is the narrow P-wave Z)+ meson with J p = 1 + . The "upper- 
vertex" production of D~£[ from W + — *> cs hadronization is shown in Figure |l](a). In addition, 
mesons can be produced from "lower-vertex" decays b — *> cud with the creation of an ss 
quark pair, as shown in Figure |l|(b). This produces right-sign D mesons; however, the decay 
rate is expected to be small. Throughout this paper charge conjugate states are implied. 

Continuum Df x production has been thoroughly studied (T|]. The is just above 
the D*K mass threshold and decays dominantly into D*°K + and D* + K°. Other possible 
decay channels are negligible: D^+tc due to isospin conservation, D^ + (nrr) due to OZI 
suppression 0, DK or Dfn° due to angular momentum and parity conservation, and .OW + 7 
due to the small radiative decay rate. 



II. DATA SAMPLE AND EVENT SELECTION 

The data used in this analysis were selected from hadronic events collected by the CLEO II 



detector at the Cornell Electron Storage Ring (CESR). The CLEO II detector |10[ is a 
large solenoidal detector with 67 tracking layers and a Csl electromagnetic calorimeter that 
provides efficient ir° reconstruction. The data consist of an integrated luminosity of 3.11 
fb _1 at the T(4S) resonance, corresponding to 3.3 x 10 6 BB events. To evaluate non-BB 
backgrounds we also collected 1.61 fb _1 of "continuum" data 60 MeV below the T(4S') 
resonance. 

The inclusive B — > D^X decay is studied by reconstructing the decay channels — > 
D*°K+ and D*+K° s using the decay modes D*° -> D°vr and D* + -> D°ir+. The D° 
is reconstructed using the decay modes D° — ► K~n + and K~7^ + tt°. Hadronic events are 



required to satisfy the ratio of Fox- Wolfram moments [11 R2 = H2/H0 < 0.3 to reduce the 
background from continuum events. 

Charged tracks, except pions from Kg decays, are required to be consistent with coming 
from the primary interaction point. Charged kaon and pion candidates are identified using 
specific ionization (dE/dx) and, when available, time-of-flight (TOF) information. For kaon 
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identification, we consider the relative probability for a charged track to be a kaon, IZk = 
Vk/CPk + Vk + V p ), where V is the \ 2 probability for a given particle hypothesis. The 
requirement on IZk depends on the decay mode of interest. Pion candidates are identified 
by requiring the dE/dx and, when available, TOF information to be within 3 standard 
deviations (a) of that expected for pions. We select K s candidates through the decay to 
7r + 7r~ by requiring a decay vertex displaced from the primary interaction point and a K s 
invariant mass within 10 MeV/c 2 of its nominal value. We reconstruct 7r° candidates through 
the decay to 77 by requiring candidates to have an invariant mass within 2.5 standard 
deviations (a ~ 5 MeV/c 2 ) of the nominal 7r° mass. 

The K~7t + and K~ir + 7r° combinations are required to have a kaon identification of IZk > 
0.5 and 0.7, respectively, and an invariant mass within 15 and 25 MeV/c 2 (~ 2a) of the 
nominal D° mass, respectively. In addition, we select regions of the D° — > K~tt + tt° Dalitz 
plot to take advantage of the known resonant substructure | 12| . For the Df x — > D*°K + 



mode, the Dalitz cut reduces the signal efficiency by 40% and the background by 80%. We 
relax the Dalitz cut for the D* + K s mode since the combinatoric background is substantially 
lower. 

The D* + — > D°7T + candidates are required to have a mass difference M(D°tt + ) — M(D°) 
within 1.5 MeV/c 2 (~ 2a) of the nominal value of 145.4 MeV/c 2 , where M(X) is the recon- 
structed invariant mass of X. Similarly, the D*° — > D°tt° candidates are required to have 
a mass difference M(D n°) — M(D°) within 1.5 MeV/c 2 (~ 2a) of the nominal value of 
142.1 MeV/c 2 . To form Df x candidates charged kaons are combined with D*° candidates 
and K s 's are combined with D* + candidates. Since the primary kaons from Dg t — > D*°K + 
decays have low momentum, we can impose a stringent IZk > 0.9 requirement on the K + 
with negligible loss of efficiency The Df x candidates are required to have a scaled momen- 
tum x„ = p D + I \/El — M 2 + < 0.45, which is the kinematic limit for B — > D^X decays. 

(We ignore the negligible contributions from b —> u decays.) Upper- vertex production 
results in a maximum x p of 0.35, and this requirement is imposed when determining the 
decay constant. The decay channels with 7r°'s in the final state often have multiple Df t 
candidates per event. We select the candidate with the highest x 2 probability of being a 
Dg X , which is derived from the invariant masses of the reconstructed ir°, D° and D* mesons. 



III. RAW YIELDS 

The D+ x signal is identified using the D*K mass difference, AM X = M(D*°K + ) — 
M(D*°) - M K + and AM 2 = M(D* + K° S ) - M(D* + ) - M K o, where M K + and M K o are 
the known masses [|l|. The D*K mass difference signal has a resolution that is two to four 
times smaller than the corresponding signal in the reconstructed D*K invariant mass dis- 
tribution. The AMi and AM2 distributions are shown in Figure |2], where the D° — ► K~tt + 
and K~tt + tt° modes have been added together. The data is fit with a Gaussian signal and 
a threshold background function. The Gaussian width is fixed to that expected from a 
GEANT-based Monte Carlo simulation |TJ[ (a = 2.4 — 3.6 MeV/c 2 , depending on the mode) 
and the mean is fixed to the measured Df x mass difference from continuum data (AMi ~ 35 
MeV/c 2 and AM 2 « 27 MeV/c 2 .) We observe 42 ± 14 signal events in the D*°K + mode and 
9 ± 6 events in the D* + K s mode. 
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However, when the D*°K + candidates are further subdivided into the D° — > K~ir + and 
K~7i + tt° decay channels there is a discrepancy in the D^ x yields. As shown in Figure ^ 
we observe 10 ± 8 signal events in the AM\ distribution for the D° — > K~tt + channel and 
33 ± 12 signal events for the D° — > K~7T + 7r° channel. After accounting for branching 
fractions and efficiencies, discussed below, this results in a 2.2a discrepancy in the D*°K + 
rates between the two D° modes. We cannot rule out the fact that background sources may 
be contributing a false signal in the D° — > K~tt + tt° channel, but not in the D° — > K~n + 
channel. However, no such mechanism has been uncovered. To be conservative, we choose 
to quote only an upper limit for the decay B — > D^X. 

Since the Df x reconstruction efficiency increases rapidly with x p and the Df x momentum 
distribution from B decays is not known, we compute the inclusive B — > D^X branching 
fraction by dividing the data into four equal regions of x p from 0.05 to 0.45 and summing 
the efficiency corrected yields. The Df x — > D*°K + and D* + K° branching fractions are equal 
according to isospin, and their ratio has been measured to be within 30% of unity |T4J. We 
measure the branching fraction B -> D^X to be (0.77 ± 0.22)% from the D*°K + mode 
and (0.28 ± 0.37)% from the D* + K$ mode, where the error is statistical only. The two 
measurements are statistically consistent. The x p distribution for our candidates is 
shown in Figure 



IV. CROSS-CHECKS 

Several cross-checks, shown in Figure [5], were performed to corroborate the validity of the 
signal. The scaled continuum background from data after satisfying all selection cuts is 
negligible, and there is no excess in the AM\ signal region (3 ±5 events). The uncertainty in 
the continuum Df x contribution is included in the systematic error. There is also no evidence 
of peaking in the AMi signal region for wrong-sign D*°K~ combinations (0 ± 9 events), D° 
mass sidebands (5 ± 5 events), and D*° mass sidebands (—4 ± 6 events). 

We have also searched for the D° signal from DJ X — > D*°K + candidates in the AM\ signal 
region, | AM\— 35 MeV/c 2 | < 10 MeV/c 2 , by relaxing the D° mass cut and histogramming the 
invariant mass of all K~ir + and K^tt + it combinations that satisfy the remaining selection 
criteria. In events with multiple candidates per D° decay mode we select the candidate 
with the highest x 2 probability, which is derived from the reconstructed ir° and Df± masses. 
We observe 100 ± 15 D° events. However, there are also real D 0, s in the random D*°K + 
combinations under the peak; after a AM\ sideband subtraction the D° invariant mass 
spectrum yields 44 ± 18 events (see Figure §(a)). This is consistent with our — ► D*°K + 
yield in Figure 0. 

Similarly, we have studied the D*° signal from Df x — > D*°K + candidates in the AM± 
signal region. We observe 59 ± 15 D° events. As in the D° case there are also real D* 0, s in 
the random D*°K + combinations under the DJ X peak. After a AM\ sideband subtraction 
the D*° mass difference spectrum yields 25 ± 18 events (See Figure |6](b)), consistent with 
our Dti -> D*°K+ yield. 

Finally, we have studied the production from continuum e + e _ — > cc events. The 
selection criteria is similar to that used to find from B decays, but since continuum 
charm production has a hard fragmentation, we require x p > 0.5. In addition, we remove 
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the i?2 < 0.3 cut, relax the charged kaon identification to TZk > 0.1, and remove the 
Dalitz cut for D° — > K^-k + -k° . The mass difference distribution for D*°K + and D* + K$ 
combinations are shown in Figure [?|, where the D° — > K~7t + and K^tt + tc modes have been 
added together. We extract the signal by fitting the data with a Gaussian signal and 
a threshold background function. The Gaussian width is fixed to the value predicted by 
Monte Carlo (2.1 MeV/c 2 ), and the mean is allowed to float. We observe 222 ± 19 events 
in the — > D*°K + mode with a mass difference of 35.0 ± 0.2 MeV/c 2 (statistical error 
only), and 101 ± 11 events in the —>■ D* + K° s mode with a mass difference of 27.5 ± 0.3 
MeV/c 2 . The results are consistent with the previous CLEO analysis W% - 



V. SYSTEMATIC ERRORS AND FINAL RESULTS 

There are several sources of systematic error. We assign a systematic error of 16% 
to account for the 2.2<r discrepancy between the D*°K + rates for the D° — > K~n + and 
K~7t + 7t modes. This accomodates different methods of computing the weighted average of 
the B — > D^X branching fraction from the four separate decay chains. Uncertainties due to 
reconstruction efficiencies include 1.5% per charged track, 5% per 7r°, 5% for slow pions from 
D* , and 5% for Kg. We also include systematic errors of 7% for Monte Carlo statistics, 5% 
for kaon identification and the Dalitz decay cut efficiency, 4% for uncertainties in the yield 
for x p < 0.05, and 8% for uncertainties in the continuum Dji contribution that passes our 
selection criteria. The total systematic error is 24%. 

Averaging the D*°K + and D* + K° s modes together, we obtain B(B — > D^X) = (0.64 ± 
0.19 ± 0.15)%. Since the signal is observed largely in only one decay mode — > 
D*°K + with D° — > K~7t + 7t°, and since there is a discrepancy between this mode and the 
corresponding mode involving D° — > K~ir + , we instead prefer to quote an upper limit on 
the branching fraction to be B < 0.95% at the 90% C.L. ]I5| This decay rate limit is small 
relative to the total rate expected for B -> D^D^KX of about (7.9 ± 2.2)% from the 
wrong-sign D meson yield in B decays 0. This is not surprising considering the cs system 
has appreciable phase space beyond the D~^ x mass ||. Also, CLEO's |16| recent observation 



of exclusive B — > D^D^*'K decays shows that the D^K invariant mass distribution lies 
mostly above the mass. 



VI. D+ DECAY CONSTANT 



si 



Measurement of the B — > D^X decay rate also provides an estimate of the decay 
constant, f D + , assuming that the comes dominantly from upper- vertex decays. The 
inclusive decay rate for B mesons into ground state or excited mesons can be calculated 



assuming factorization |17 



T(B -> D S X) 



G 2 p\V C bV C s\ 
16tt 



Mla\f Ds I{x,y) 



where a\ is the BSW [|TB[ parameter for the effective charged current, and I(x,y) is a kine- 
matic factor with x = M|, /M 2 and y = M^/M 2 . For scalar or pseudoscalar D s mesons, 
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I(x,y) — J(l — x — y) 2 — 4xy(l — x — 2y — xy + y 2 ), and for vector or axial- vector D s mesons, 



I(x, y) = y (1 — x — y) 2 — Axy{l + x — 2x 2 — 2y + xy + y 2 ). 

We have tightened the x p requirement to x p < 0.35 since this is the kinematic limit 
for upper-vertex B — > D^DX decays. The production of ground state and excited D + 
mesons from lower-vertex decays such as B — > D^KX is expected to be suppressed. This is 
certainly true for B — > D + X decays where the fraction of D + produced at the lower-vertex is 
measured to be 0.172 ± 0.079 ± 0.026 [HJ. Moreover, there is no evidence of D + ± production 
in the region x p = 0.35 — 0.45 where lower-vertex production is likely to occur (see Figure^.) 

With the assumption f D + = f D *+ we can extract f D + from the ratio of inclusive rates, 

B(B - D+X) _ T(B - D+X) ^ Q ^ ( f D + \ 2 



B(B — > D+X) T(B -»■ D+X) + F(B -»■ D*+X) ' V/_ 



Many systematic errors cancel in the ratio. When computing the decay constant from the 
above equation, we use (75 ±25)% of the measured B — > D^X branching fraction to account 
for uncertainties in the upper and lower vertex contributions to D~£[. This accomodates the 
excess of B — > D^X candidates observed at low x p < 0.15 as seen in Figure £|. From our 
upper limit on B ^ D^X and CLEO's measurement of B(B -> D+X) = (12.11 ± 
0.39 ± 0.88 ± 1.38)%, we derive f D + I f D + < 0.40 at the 90% C.L. The central value is 

si s 

f D + I f D + = 0.29 ± 0.06 ± 0.06, where the first error is due to the total error in the inclusive 

si s 

B — > D+X and B — > D^X branching fractions, and the second is the uncertainty in the 
non-factorizable and lower-vertex contributions to the B — > DJyX decay rate. Using the 



measured value of f D + = 280 ±40 MeV p0| gives f D + = 81 ±26 MeV which corresponds to 



an upper limit of f D + < 114 MeV. This limit accomodates the prediction of f D + = 87 ± 19 



MeV by Veseli and Dunietz [21] 



VII. CONCLUSIONS 

In summary, we have searched for B mesons decaying into the P-wave .0^(2536) meson. 
The upper limit of B(B — > D+ X X) < 0.95% at the 90% C.L. accounts for at most only a 
fraction of the total wrong-sign B — > DX rate. Assuming factorization, the decay constant 
f D + is at least a factor of 2.5 times smaller than the decay constant for the pseudoscalar 

si 

Dt- 
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FIG. 1. Feynman diagrams for (a) B — ► D^X decays producing at the upper-vertex and 
(b) B — ► decays producing D~ x at the lower-vertex. 
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FIG. 4. The efficiency corrected yield for our B — > D^X candidates as a function of the Df t 
scaled momentum x p . The kinematic limit from upper-vertex and lower-vertex B — > D^[X decays 



is x p < 0.35 and x p < 0.45, respectively. 
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FIG. 5. The normalized D*°K + mass difference distributions from (a) continuum events, (b) 
D*°K~ "wrong-sign" combinations, (c) D° mass sidebands, and (d) D*° mass sidebands. 
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FIG. 6. (a) The invariant mass distribution for K~ir + and K^tt + tt° combinations from D*°K + 
candidates in the AM\ signal region, after sideband subtraction, (b) The D*° mass difference 
distribution from D*°K + candidates in the AMi signal region, after sideband subtraction. 
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FIG. 7. The mass difference distribution for (a) D*°K + and (b) D* + K$ candidates from con- 
tinuum e + e _ — > cc events. 
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